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Abstract 

In the last few years considerable observational resources have been devoted to 
study the thermal emission from isolated neutron stars. Detailed XMM and Chandra 
observations revealed a number of features in the X-ray pulse profile, like asymmetry, 
energy dependence, and possible evolution of the pulse profile over a time scale of 
months or years. Here we show that these characteristics may be explained by a 
patchy surface temperature distribution, which is expected if the magnetic field has 
a complex structure in which higher order multipoles contribute together with the 
dipole. We reconsider these effects from a theoretical point of view, and discuss their 
implications to the observational properties of thermally emitting neutron stars. 
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1 Introduction 



The seven X-ray dim isolated neutron stars (XDINSs) discovered so far (see 
e.g. Treves et al. 2000, Motch 2001 for a review) offer an unprecedented oppor- 
tunity to confront present models of neutron star (NS) surface emission with 
observations. These objects play a key role in compact objects astrophysics 
being the only sources in which we can have a clean view of the compact 
star surface. In particular, when pulsations and/or long term variations are 
detected we can study the shape and evolution of the pulse profile of the ther- 
mal emission and obtain information about the thermal and magnetic map of 
the star surface. So far, X-ray pulsations have been detected in four XDINSs, 
with periods in the range 3-11 s. In each of the four cases the pulsed fraction 
is relatively large (~ 12% — 35%) and, counter-intuitively, the softness ratio 
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is maximum at the pulse maximum (Cropper et al. 2001, Haberl et al. 2003). 
Spectral lines have been detected in the soft X-rays and the line parameters 
may change with spin pulse. In addition, often the X-ray light curves appears 
to be slightly asymmetric and a gradual, long term evolution in both the X- 
ray spectrum and the pulse profile of the second most luminous source, RX 
J0720. 4-3125 has been recently discovered (De Vries et al. 2004). 
All these new findings represent a challenge for conventional atmospheric mod- 
els: the properties of the observed pulse profiles (large pulsed fraction, skew- 
ness, and possibly time variations) cannot be explained by assuming that the 
thermal emission originates at the NS surface if the thermal distribution is 
induced by a simple core-centered dipolar magnetic field. On the other hand, 
it has been realized long ago that two effects can contribute at least to achieve 
relatively large pulsed fraction (up to 20%) : 1) radiative beaming (Pavlov et 
al. 1994) and 2) the presence of quadrupolar components in the magnetic field 
(Page and Sarmiento 1996). Here we present magnetized atmospheric mod- 
els computed assuming a quadrupolar magnetic field geometry and show how 
they can account for some of the general characteristics of the observed X-ray 
lightcurves (see Zane and Turolla 2004, in preparation for further details). 



2 Getting to grips with the neutron star surface 

2.1 Computing the light curve 

In order to compute the phase-dependent spectrum emitted by a cooling NS 
as seen by a distant observer, we basically perform three steps. First, we as- 
sume that the star magnetic field possesses a core-centered dipole+quadrupole 
topology, B = B dip + B quad , where B quad = Y%=o & B JLi- The polar compo- 
nents of ~B dip and of the five generating vectors B^ ad are reported in Page and 
Sarmiento (1996). The NS surface temperature distribution is then computed 
using the simple expression T s = T p \ cos a] 1 / 2 , where T p is the polar tempera- 
ture and a is the angle between the field and the radial direction, cos a = B • n. 
Second, we compute the local spectrum emitted by each patch of the star sur- 
face by using fully ionized, magnetized hydrogen atmosphere models. 1 Besides 
surface gravity, this depends on both the surface temperature T s and magnetic 
field, either its strength and orientation with respect to the local normal. We 
introduce a (6, 0) mesh which divides the star surface into a given number of 
patches. The atmospheric structure and radiative transfer are then computed 
locally by approximating each atmospheric patch with a plane parallel slab. 

1 We caveat that partial ionization effects are not included in our work. Bound 
atoms and molecules can affect the results, changing the radiation properties at 
relatively low T and high B (Potekhin et al. 2004). 
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Third, we collect the contributions of surface elements which are "in view" 
at different rotation phases (see Pechenick et al. 1983, Lloyd et al. 2003). We 
take the neutron star to be spherical (mass M, radius R) and rotating with 
constant angular velocity uj = 2n/P, where P is the period. Since XDINs are 
slow rotators (P 10 s), we can describe the space-time outside the NS in 
terms of the Schwarzschild geometry. Under these assumptions, for a fixed 
polar temperature, dipolar field strength and surface gravity M/R, the com- 
puted light curve depends on seven parameters: 6, = qi/B dip (i = 0, ...4), the 
angle x between line of sight and spin axis, and the angle £ between magnetic 
dipole and spin axis. 



2.2 Studying light curves as a population and fitting the observed pulse shapes 

Given this multidimensional dependence, an obvious question is whether or 
not we can identify some possible combinations of the independent param- 
eters that are associated to particular features observed in the pulse shape. 
Particularly promising for a quantitative classification is a tool called princi- 
pal components analysis (PCA), which is concerned with finding the minimum 
number of linearly independent variables z p (called the principal components, 
PCs) needed to recreate a data set. In order to address this issue, we divided 
the phase interval (0 < 7 < 1) in 32 equally spaced bins, and we computed a 
population of 78000 light curves by varying x, £ an d h (i = 0, ..4). By applying 
a PCA, we found that each light curve can be reproduced by using only the 
first ~ 20 more significant PCs and that the first four (three) PCs account for 
85% (72%) of the total variance. However, due to the strong non-linearity, we 
find so far difficult to relate the PCs to physical variables. 
Nevertheless, the PCA can be regarded as an useful method to provide a link 
between the various light curves, since models "close" to each other in the PCs 
space have similar characteristics. From the PCA we obtain the matrix CV, 
such that Zi = C^yj, where yj is the observed X-ray intensity at phase 7^ and 
Zi is the i-th PC. Therefore, we can compute the PCs corresponding to each 
observed light curve and search the models population for the nearest solution 
in the PCs space (see fig. 1, left). This in turn provides us a good trial solution, 
which can be used as a starting point for a numerical fit. The quadrupolar 
components and viewing angles are treated as free parameters while the polar 
values of T p and B>di P are fixed and must fall in the domain spanned by our 
atmosphere model archive. 2 Our preliminary results are illustrated in Fig. 1, 
2, and 3, and refer to B dip = 6 x 10 12 G, logT p (K) = 6.1 - 6.2. 



2 In general this will not contains the exact values inferred from spectral observa- 
tions of XDINSs. However, we have checked that a fit (albeit with different values 
of the quadrupolar field) is possible for different combinations of B^ip and T p in the 
range of interest. 
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3 Summary of results 



As we can see from Fig. 1, 2 and 3, the broad characteristics of all the XDINS 
light curves observed so far are reproduced when allowing for a combination 
of quadrupolar magnetic field components and viewing angles. However, al- 
though in all fit exists, we find that in general it may be not unique. 
This model has not a "predictive" capacity in providing the exact values of 
the magnetic field components and viewing angles: this is why we do not fit 
for all the parameters in a proper sense and we do not derive parameter un- 
certainties or confidence levels. The goal is instead to show that there exist 
at least one (and more probably more) combination of parameters that can 
explain the observed pulse shapes, while this is not possible assuming a pure 
dipole configuration. In the case of RX J0720. 4-3125, preliminary results show 
that the pulse variation observed between rev. 78 and rev. 711 can not be ex- 
plained by a change in viewing angle only (as it should be if NS precession 
is invoked, de Vries et al. 2004) or magnetic field only. Instead, a change in 
all quantities (quadrupolar components and viewing angles) is needed. 3 Aim 
of our future work is to reduce the degeneracy either by performing a more 
detailed statistical analysis on the models population and by refining the best 
fits solutions by using information from the light curves observed in different 
color bands and/or from the spectral line variations with spin pulse. 
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• Fig. 1. Left: The computed population of 78000 light curves plotted against 
the first 3 principal components. Axis are arbitrary; yellow symbols mark 
the position of the EPIC-PN light curves of the XDINSs (see next captions 
for references; both rev. 78 and rev. 711 are shown for RX J0720. 4-3125), 
red symbols mark the light curves computed assuming a pure dipolar field. 
The PCs representation (limited to the first three PCs, which alone account 
for the 72% of the total variance) of the observed light curves falls within the 
domain spanned by the quadrupolar model representations; this is why at 
least one fitting solution can be found. Right: Fit of the EPIC-PN (0.12-0.7 
keV) smoothed light curve of RX J0420.0-5022 (Haberl et al. 2004, rev. 570). 
Dashed line: trial solution as inferred from the closest model in the PCs 
space; solid line: best fit solution. The best fit parameters are: b = —0.48, 
h = 0.02, b 2 = -0.25, 6 3 = 0.35, 64 = -0.20, f = 39.9°, x = 91.2°. 

• Fig. 2. Same as in the right panel of Fig. 1 for two further INSs. Left: Fit 
of the EPIC-PN (0.12-1.2 keV) light curve of RX J0806.4-4123 (rev. 618, 
Haberl et al. 2004). Here: b = 0.39, b x = -0.37, b 2 = 0.12, 6 3 = -0.13, 
b 4 = 0.49, f = 0.0°, x = 59.2°. Right: fit of the EPIC-PN (0.12-0.5 keV) 
light curve of RBS 1223 (rev. 561, Haberl et al. 2003). Here: b = 0.21, 
h = -0.02, b 2 = 0.59, b 3 = 0.53, 64 = 0.50, f = 0.0°, x = 95.1°. 

• Fig. 3. Same as in the right panel of Fig. 1 for the fit of two EPIC-PN (0.12- 
1.2 keV) light curves of RX J0720.4-3125 (both from De Vries et al. 2004). 
Left (rev. 78): b = 0.36, b x = 0.43, b 2 = -0.16, 63 = -0.16, b 4 = -0.39, 
f = 0.0°, x = 68.1°. Right (rev. 711): b = 0.45, b x = 0.49, b 2 = -0.06, 
63 = -0.08, h = -0.26, f = 0.0°, x = 87.6°. 
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Fig. 1. 




Fig. 2. 




Fig. 3. 
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